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We calculate the rapidity dependence of the transverse double-spin asymmetry for the Drell-Yan 

0^ ■ process to next-to-leading order in the strong coupling. Input transversity distributions are obtained 

0^ ' by saturating the Soffer inequality at a low hadronic mass scale. Results for the polarized Bnl- 

O^ , Rhic proton-proton collider and the proposed Hera-N fixed-target experiment are presented, and 

the influence of the limited muon acceptance of the detectors on measurements of the asymmetry is 

P3 ' studied in detail. 

^' 

p. 1 One of the major goals of the forthcoming spin programme at the Rhic polarized proton-proton collider [n| is a first 

measurement of the twist-2 transversity distribution (5q(x,/i^) [g|, which is theoretically as important as the well- 

fS| known unpolarized and longitudinally polarized parton densities q{x,fp) and Ag(a;,/i^), respectively. An important, 

^ non-trivial model-independent restriction on the size of Sq{x, fj?) derives from Soffer's inequality 0], which states that 

o . 

^ : \5q{x, /i^)! < i [q{x, /i^) -I- Ag(x, m')] , (1) 

(n: 

^^ ■ and similarly for antiquarks. It was shown to be preserved by next-to-leading order (NLO) DGLAP evolution in 
"reasonable" factorization schemes, among them the MS-scheme [4-6]. 
In a previous NLO analysis, see Q for details, we have derived an upper bound [by saturating (^] for the total 

1-^ ' transverse double-spin asymmetry AttIM) for Drell-Yan dimuon production of mass M . It turned out, however, that 
t^. Att{M) is not very sensitive to the shape of 5q. In addition the angular acceptance of the detectors was assumed to 
O , ' be constant, i.e., independent of the dimuon rapidity y, which can only be a rather crude approximation to the real 
1^ , experimental conditions. Therefore, in order to better suit the experimental needs, we extend the analysis of Iq] in 

T^ • this note and study the y dependence of Att- 
►>. The transversely polarized Drell-Yan cross section, A5a = (du^^ — dcr^-'') /2, is given as a double convolution of 

• ^ 1 transversity distributions with the corresponding transversely polarized partonic cross section: 

X; 

■ ^- '• d]i^ = ^^~'//^^//^^ [Sq{x,,^l)Sq{x,,^l)+6qix,,^,%)Sq{x„^l)] ^^^ , (2) 

I^F being the factorization scale. The effective charge Cg also contains the electroweak effects from Z" exchange and 
7^° interference; see, e.g., Eq. (20) of [p|. i/ denotes the rapidity of the dimuon pair, and cj) is the azimuthal angle of 
one muon, with (j) = in the direction of positive transverse spin of the incoming protons. The variables x^, a;^ in 
(0) are related to y and the Drell-Yan scaling variable t = A'P / S by x\ = y/rey and a;^ = \pTe~^ , so that the region 
y > (y < 0) is mainly sensitive to x\ [x^. To lowest order (LO) x\ and x\ coincide with the momentum fractions 
carried by the incident partons. Indeed, one has at LOQ 

^^'"'^ ^"' cos{2<f>)&[x^-x\)&[x2-x\). (3) 
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dA'fdyd^ 95M 
In the MS-scheme, the NLO [©(as)] correction to (||) reads 



^The formula $(0) = 1 below Eq. (15) in |i] should read ^{<f>) = 2. 



d5a(i)^^s 



2a a,{^i%) 



4:T{xiX2 + t) 



dMdydcj) 9SM'^"' 2tt xiX2{xi + x°i){x2 + x^) 
X \s{xi-x1)S{x2-x°) 

1 



(K cos(20) 



1 2 (l-2:?)(l-2;") 7r2 
4 T 4 



+5{xi -a;?) 



(X2 - 0:^)+ 



In 



2a;2(l-x?) , /ln(a;2-x^) 



-(xa+xO) 



1 



■In- 



X2 - x'^ J ^ X2- a;^ 



X2 



+ 



(Xi ^X\){X2+X\) 



3 In 



zia;2+r 



2[{X,-X°,){X2-X^,)], 



In 






2 r 



S{xi — x")(5(x2 — a;2) 



{xix'^+X2X^y {Xi - X'(){X2 - x'^) 

4 



2 



(5(xi 






(X2 - X§)h 



where /i^^ is the renormahzation scale (for simpHcity we always take ^r = i-ip = M) and (i = 1, 2) 



dxif(xi) 



1 



{xi - a;°)+ 



\^ f{x,)~J{x^l) 



+ [l--2], (4) 



(5) 



Equation (Q) is obtained by a suitable factorization-scheme transformation of the corresponding result of [Q , which 
was calculated taking the gluon off-shell in the process qq — *■ fi^ pT g in order to regularize its coUinear divergences. 
The corresponding results for the unpolarized NLO cross section dcr = (dcr^^ -f dcr^^) /2 can be found in gQ 

In order to increase the observable rates, we will integrate the unpolarized cross section over (^, whereas in the 
polarized case we add each quadrant with a different sign. Thus, the rapidity dependent asymmetry will be defined 
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where Mq^i denote the limits of some suitable bin in invariant mass. Following closely our previous study g on the 
total (i.e. y-integrated) Drell-Yan cross section, we will try to estimate upper hounds on Att by assuming that the 
equality in (|^) holds[| at a low hadronic mass scale jjlq ~ 0(0.6 GeV), see Q for more details. We should emphasize 
that the sign of the asymmetry cannot be predicted in this way, because only the absolute value of 5q enters Soffer's 
inequality. This also means that all possible combinations of signs in Eq. (|^) must be checked so as to obtain the 
maximal absolute value of Att- In our case, choosing all signs to be the same always yielded the largest results. 

Figures ^ and || show the "maximally possible" dda/dy and Att in LO and NLO for ^/S = 500 GeV at Rhic and for 
Ehcum = 820 GeV, corresponding to v^ = 39.2 GeV, at Hera-N, respectively. We have integrated over M in (||) as 
indicated in the figures, avoiding masses smaller than 5 (4) GeV for Rhic (Hera-N), where a large background from 
charmed-meson decay is expected. Very similar results as in Fig. |l| are obtained for y/S = 200 GeV and C = 320 pb^^ 
at Rhic when restricting M to be in the range 5 — 9 GeV. The QCD corrections to the polarized cross section turn 
out to be largest in the fixed-target regime, whereas the asymmetry receives the largest corrections at higher energies. 
In most cases the NLO contributions are sizeable and should be included for a meaningful comparison with future 
data. We note in passing that we found that the dependence of the results on ^ji and fip is greatly reduced at NLO. 

In Figs. H and || we also display the statistical errors expected for such measurements of Att- Here, we try to 
estimate the influence of detector cuts on the error, which could be rather crucial for making realistic predictions. For 
instance, if the muon detectors have limited angular coverage, one or both of the muons might escape detection just 
for geometrical reasons, and the event is lost. In the case of the Rhic detector Phenixq, the endcaps will be able 



The result of M is not given in the conventional MS-scheme; however the translation can be easily made. 

^In H we actually did not saturate the total quark distributions, but only their valence component at the input scale /iq. As 
was pointed out in M, this is, strictly speaking, not the statement of the Softer inequality. A careful numerical check however 
reveals that none of our results in M is altered if one saturates the full quark distributions instead of the valence ones. 

*We only calculate acceptance corrections for Phenix, since the other major Rhic detector. Star, cannot detect muons, but 
only electrons. Electron pair production does not seem as promising as muon pair production, as a very detailed study of the 
background is required in that case. 



to identify muons with 1.2 < |2/^±| < 2.4; an additional cut on the muon momentum, |A:| > 2 GeV, will probably be 
necessary to get rid of unwanted background. Central rapidity muon detector arms, which would cover |y^±| < 0.35 
(even though for only half of the azimuth) , were proposed but will not be realized [Q . Nevertheless, we have also 
studied the impact that they would have had on the achievable experimental accuracy. In order to calculate the 
relevant acceptances, the momenta of the outgoing muons must be known. However, they cannot be reconstructed 
from the kinematic variables Af , y, and cj) introduced above, since M and y refer to the dimuon system and (j) is only 
one of the angles describing the direction of one muon. Therefore, one has to consider a more differential cross section, 
like 






TP 



where ({M, fc^, 0) — kx {2 — Ak^/M^), 5C,{M, kr, 4>) — 4cos(20)fc^/M^, and kr is the transverse momentum of the 
muons. The LO acceptance curve for the measurement of, say, they-dependence of the cross section or the asymmetry 
Atti can then be obtained by dividing the results based on Eq. ([7|), after implementation of appropriate cuts on t/^±, 
by the full LO result, i.e., the one integrated over all kx and already used in Figs. |l] and ||. Of course one could 
have also extended the acceptance analysis to NLO, where the muons are no longer back-to-back and the possibility 
arises that both muons go into the same hemisphere of the detector. However, we believe that a LO estimate for the 
acceptance is good enough to get a rough quantitative understanding of the influence of limited detector coverage on 
the statistical error. 

Figure H shows the acceptances for muon identification in the endcaps only and for the endcaps plus central 
detector arms. Note that the unpolarized acceptances e differ from the polarized ones Se as a result of the different 
fcT-dependences of the corresponding cross sections (0). The results for VS* = 200 GeV and 500 GeV turn out to be 
almost the same, because we used the same lower limit for the dimuon mass M in both cases. According to Fig. 0, 
the acceptance for the central rapidity region y sa 0, where each endcap or each central arm detects one muon, is 
considerably smaller than for the large rapidity region, where both muons hit the same endcap. Also, the ratio of 
"polarized-to-unpolarized acceptance" is smaller than unity in the former case and larger than unity for the latter. 
This means that the experimentally measured asymmetry will be smaller at y « 0, but somewhat enhanced at large 
y as compared to the values given in Fig. n^. We also see that the addition of muon identification in the central arms 
would yield a much larger acceptance at small and intermediate dimuon rapidities than found for the "endcaps only" 
scenario. 

At the moment, Hera-N only has the status of a fairly general proposal for a fixed-target pp spin experiment at 
Hera [ pT| . Thus, nothing specific is known yet about appropriate kinematical cuts. In our analysis we try to use 
reasonable values for the kinematical coverage, keeping in mind that the true detector could look significantly different 
in case it will ever be built. We use ±700 mrad for the horizontal and ±160 mrad for the vertical opening angle, while 
the beam pipe is assumed to cover ±10 mrad. Such a detector would have much larger acceptances than Phenix, as 
can also be seen in Fig. y. 

Exploiting our LO estimates of the acceptances e and (5e, we are now in a position to calculate the expected 
statistical errors on the asymmetry. Here we assume that it makes sense to adopt our LO acceptances curve also for 
the NLO calculation; see our discussion above. The statistical error of the "measured" asymmetry, i.e., a/ter correction 

for acceptance, is then just given by l/V^-JC J eda where V denotes the degree of polarization of each beam, C is 
the integrated luminosity, and the integration goes over the bin under consideration. In order to consistently match 

the error bars to Figs. || and ||, we obviously have to weigh them by the ratio "K- — / '^ . 

The statistical errors show the same features for both Rhic energies. A measurement in the central rapidity region 
will hardly be possible, even if the central muon detector arms are added. Statistical errors at large rapidities do 
not depend on the presence of central rapidity muon detection (see Fig. y), and prospects look slightly better here. 
The larger rates for ^/S — 500 GeV are compensated by a smaller asymmetry so that, for both ^/S — 200 GeV and 
^/S = 500 GcV, the relative statistical error is about 40% at large y. Note that we also include the events with 
negative rapidity for the calculation of the error bars, since the results are symmetric in y. The situation for Hera-N 
is somewhat better, with relative errors of about 30%, and more possible bins. This is mainly due to the much larger 
asymmetry in the fixed-target regime. However, for all this we should keep in mind that the asymmetries we show 
have been obtained assuming saturated 6q^s at a low scale. If the saturation were only at, say, the 50% level, then all 
asymmetries would have to be scaled down by a factor 4, and no measurement would be possible. 

Clearly, the restriction in angular acceptance expressed by Fig. will also leave its footprint for the y-integrated, i.e., 
the total, Drell-Yan cross section. In other words, we have to reinspect our predictions made in H for this quantity, to 
see whether there is any dramatic change concerning the statistical accuracy of a possible measurement of Att{M). 



On the left-hand sides of Figs. H and || we show the unpolarized and polarized acceptances for the total dimuon cross 
section for the Rhic energy of yS — 200 GeV and the Hera-N situation. In the case of Rhic, we distinguish again 
between the "endcaps only" and the "endcaps plus arms" options. The general trend is that the acceptances are 
rather low for Rhic (Phenix) and decrease with increasing M after reaching a peak at a quite low M-value. Under 
our assumed conditions for Hera-N, the acceptance turns out to be much higher and fairly independent of M. On the 
right-hand sides of Figs. |4| and ^ we redisplay our findings for Att{M) of Figs. 3 and 4 of Rcf. Q, but now with the 
more realistic error bars based on our considerations concerning the acceptance. One finds that at not too large M, 
a measurement of a non- vanishing asymmetry for the total Drell-Yan cross section still looks possible also for Rhic, 
provided the "true" transversity densities are anywhere near the ones we have modeled. Measurements at large M 
appear hopeless. The situation for ^/S = 500 GeV at Rhic is qualitatively very similar and hence not shown. Again, 
as in the case of AxTiy), Hera-N looks in a somewhat better shape. 

In conclusion, we have studied the "maximally possible" Att, resulting from saturation of Soffer's inequality at 
a low hadronic scale. It turns out that the limited muon acceptance for the Rhic experiments threatens to make 
a measurement of transversity elusive. In particular, it will be difficult, if not impossible, to measure the rapidity 
dependence of A^t, which in principle would be expected to be sensitive to the shape of 5q. At best, one data point 
at large y can be obtained, but with a large relative error. The limitation in the muon acceptance also affects the 
y-integrated cross section, so that the resulting Att{M) will also receive a substantial relative statistical error. An 
upgrade of the Phenix detector towards muon identification also in the central arms would not improve the situation 
significantly. Lower energies, in combination with better muon acceptance, seem more favorable. 
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FIGURE CAPTIONS 

Fig. 1 "Maximal" polarized cross section and asymmetry as functions of dimuon rapidity y for Rhic at vS — 
500 GeV. The error bars have been calculated for C = 800 pb^ , 70% polarization of both beams, and include 
acceptance corrections (see text). The point at low rapidity can only be obtained if Phenix is endowed with 
central muon detector arms. 



Fig. 2 Same as Fig. 1, but for Hera-N with E], 



820 GeV and C = 240 pb"!. 



Fig. 3 Acceptance curves for the detection of dimuons with the Phenix and Hera-N detectors, as functions of the 
dimuon rapidity y. The Phenix acceptances for vS = 500 GeV and Af = 5-20 GeV differ only very shghtly 
from the results shown here for the case ^/S = 200 GeV. 

Fig. 4 Dependence of the acceptances and the NLO asymmetry Att on the dimuon invariant mass, integrated over 
rapidity, for ^/S — 200 GeV at Rhic. The error bars on the right-hand side include the acceptance corrections 
and are based on £ = 320 pb~^ and V = 0.7. The outer error bars correspond to the "endcaps only" option, 
while the inner ones have been obtained assuming additional central detector arms. 

Fig. 5 Same as Fig. 0, but for \/S = 39.2 GeV, corresponding to Hera-N. 
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